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We investigate the ID Anderson-Hubbard model at half filling with box-disorder. The ground 
state phase diagram is obtained by means of real-space dynamical mean-field theory (R-DMFT) 
and the density matrix renormalization group (DMRG). We find Mott insulating and Anderson 
localized regimes as well as a strong indication of a delocalized phase for intermediate interaction and 
disorder strength within accessible system sizes. These phases are characterized and distinguished by 
qualitatively different scaling behavior of the local density of states, the energy gap in the excitation 
spectrum and the inverse participation number. 
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I. INTRODUCTION 

Disorder is omnipresent in solids due to defects and 
impurities. In his breakthrough theoretical investigation 
of non-interacting disordered systems Anderson showed 
that coherent backscattering processes due to impurities 
are responsible for spatial localization of electronic wave 
functions i The scattering probability increases with the 
impurity density and beyond a critical value leads to spa- 
tial localization of the electrons. If the states at the Fermi 
level become localized, the system becomes an Anderson 
insulator. Whereas in 3D systems the disorder strength 
has to overcome a certain threshold to drive the system 
from a metal to an insulator, in 2D and ID systems any 
nonzero disorder strength is sufficient to cause localiza- 
tion in the thermodynamic limiti 2 -^ 

On the other hand, a Mott-Hubbard metal-insulator 
transition in a pure system can be induced by electron- 
electron interaction for integer filling 4£ In contrast to 
scattering by impurities, this many-body effect induces 
correlations and tends to distribute particles uniformly. 
Numerical and analytical investigations via dynamical 
mean-field theory, exact diagonalization, quantum Monte 
Carlo and analytical calculations revealed a critical inter- 
action strength in 3D and 2D at which the phase transi- 
tion takes placed— In one spatial dimension in the ther- 
modynamic limit the Mott phase is present at any inter- 
action strengthj 3 ' 10 ' 11 

The physics resulting from the simultaneous presence 
of both effects is not yet fully understood. Although in- 
teraction and disorder both lead to metal-insulator tran- 
sitions, their action is competing: the repulsive electron- 
electron interaction favors uniform distribution of the 
particles, while disorder localizes the electronic wave 
function to a few lattice sites. 

Several investigations have been performed on the 
competition between interaction and disorder in the 
Anderson-Hubbard model: such as a perturbative renor- 



malization group (RG)^ 2 - analysis and non-perturbative 
treatments via density matrix renormalization group 
(DMRG) and dynamical mean-field theory (DMFT) 
techniques, to name just a few^ - — In these studies the 
interplay was found to affect the system in a subtle way. 
For a semi-elliptic density of states in high dimensions 
the presence of both disorder and interaction was found 
to lead to a disordered metal surrounded by an insulat- 
ing phase, the latter consisting of Anderson and Mott 
insulators continuously connected with each other i 15 ' 17 
Since in the first studies of non-interacting, weakly inter- 
acting, as well as strongly interacting 2D system s 2 ' 21 ' 22 
only insulating regions were found, it was believed for 
a long time that a metallic phase cannot occur. How- 
ever, theoretical studies via R G 13 ' 23 involving N flavors 
of electrons predicted the existence of a quantum critical 
point, at which a metal- insulator transition takes place 
in 2D. This prediction was confirmed experimentally a 
few years later i24 In ID bosonic systems a Mott insula- 
tor, a Bose glass and a superfluid phase were found ^r— 
of which the latter is analogous to the metallic phase in 
fermionic systems. However, for repulsively interacting 
fcrmions in ID perturbative RG calculations predicted 
a random antiferromagnet, where the fermions localize 
individually around randomly distributed sites ^ Never- 
theless, perturbative treatments of disorder can not be 
expected to capture the physics in the full range of cou- 
plings. In ID up to now only next-neighbor interacting 
disordered spinless fermions have been studied exactly 
via DMRG^i where, based on the phase sensitivity, a 
delocalized phase for intermediate attractive interaction 
and disorder strength was found. The ground state phase 
diagram of spinful fcrmions with on-site interactions has 
not been determined yet. 

The careful experimental analysis of these phenom- 
ena requires tunability of the disorder and interaction 
strengths. During the last decade it has become possi- 
ble to simulate theoretical models of solid state systems 
by ultracold atoms in optical lattices^ - — Due to the 
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precise control over the system parameters and the pos- 
sibility to artificially introduce different kinds of disorder, 
e.g. binar y 34 ! 35 , bichromatio^r— and speckl o 39 ' 40 , these 
systems are highly suited for detailed investigations of 
disorder phenomena. In particular, the spatial dimen- 
sionality can be easily adjusted and the localization of 
bosonic matter waves in ID 39 ' 41 and of fermions in 3D^ 
has already been observed. 

A non-perturbative investigation of a spin-i ID 
fermionic system with box disorder is the aim of this pa- 
per. By varying the disorder and interaction strengths, 
Anderson- and Mott-insulating regimes as well as der- 
ealization are found within our simulated system sizes, 
based on the real-space extension of DMFT (R-DMFT) 
and on DMRG. For their characterization, the following 
physical observables were calculated: a) the geometric 
average of the local density of states, representing its 
typical value, and its scaling behavior with the system 
size, b) the charge gap in the thermodynamic limit, c) 
the inverse participation ratio and its dependence on the 
system size. 

The paper is organized as follows: in Sec.|H]we intro- 
duce the Anderson-Hubbard Hamiltonian. The comple- 
mentary DMRG and R-DMFT methods applied in this 
paper are briefly explained in Sec. IIIII The resulting 
physical observables and a detailed analysis of the emerg- 
ing delocalized, Anderson and Mott- Hubbard regimes are 
discussed in Sec.[lV] In Sec.|V]wc conclude with a sum- 
mary of our results. 



III. METHODS 

Commonly used approaches to solve interacting quan- 
tum problems in one spatial dimension are the pertur- 
bative renormalization group (RG), the density matrix 
renormalization group (DMRG) and quantum Monte- 
Carlo (QMC). The RG is able to capture localization and 
derealization effects for repulsive as well as for attrac- 
tive interactions. However, it describes the system accu- 
rately only for small disorder strength i 12 i 43 In this paper 
we make use of DMRG which tackles strong disorder and 
strong interactions simultaneously, and allows for a de- 
termination of the phase diagram in a broad parameter 
range, similar to QMC. An alternative non-perturbative 
method for solving correlated fermionic problems is the 
dynamical mean- field theory (DMFT). In combination 
with the geometric disorder average of the local density 
of states, this method allows to detect Anderson local- 
ization as welL 15 ' 44 " — This latter approach is commonly 
referred to as typical medium theory. Within DMFT 
the self-energy is approximated to be local, which is only 
exact in infinite dimensions and is known to lead to qual- 
itatively accurate results in 3D. Its real-space extension, 
on the other hand, treats the single-particle problem - in- 
cluding disorder - exactly in any spatial dimension, while 
the self-energy becomes site-dependent 42r~— For this rea- 
son R-DMFT is expected to be superior to single-site 
DMFT in low dimensions. Here we present a quantita- 
tive analysis by a comparison of R-DMFT and DMRG 
results in ID. 



II. MODEL 

A. Hamiltonian 

Strongly correlated disordered fermions on a lattice are 
described by the Anderson-Hubbard Hamiltonian 

U = -<X)( c U Ci + 1 '< T + 4+1,0-0 + U Y1 n n n ii 

i,<7 i 

+ ^2( £ i - V) n io > (!) 
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where a 6 {t>i} labels spin, while c\ a , Ci„ and rii t(r 
are the creation, annihilation and particle number op- 
erators for an electron on site i with spin a. [/is the 
on-site interaction and t is the nearest-neighbor hop- 
ping matrix element. The on-site energies are ran- 
dom variables, each distributed independently accord- 
ing to V(si) = Q{D - |ej|) • 2/D. Here 8 is a Heavi- 
side function and D is the disorder strength. We con- 
sider a ID bipartite lattice with commensurate filling 
= (nj-j- + n.ii) = !• I 11 the homogeneous case, i.e. 
for D = 0, this Hamiltonian can be solved exactly by 
means of the Bethe ansatz^ However, advanced numer- 
ical methods are required when the on-site energies are 
random. 



A. DMRG 

Our implementation of DMRG^ 1 - - — is based on a ma- 
trix product state variational formulation ! 54 ' 55 The al- 
gorithm is formulated in the canonical ensemble (fixed 
number of particles N) and at zero temperature. It 
makes use of the SU{2) spin symmetry of the Hamilto- 
nian to significantly reduce the computational effort. The 
states are multiplet representations in the group theoret- 
ical sense, such that the actual number of different states 
is much larger than in a calculation with Abelian sym- 
metries only. We use system sizes of up to 128 sites and 
keep up to 600 states for the ground state calculations. 

When strong disorder is present, the algorithm has 
the tendency of getting trapped in excited (metastable) 
states with a slightly higher energy than the ground 
state. To solve this problem we apply nonlocal changes 
to the wave function by adding and removing a delocal- 
ized particle in a way that minimizes the energy. In most 
cases, this overcomes the barrier between the metastable 
states and the ground state such that subsequent DMRG 
sweeps will converge to the true ground state. We vali- 
date this convergence by starting DMRG from different 
initial states chosen from either the exactly solvable case 
with no interaction, the homogeneous case with no dis- 
order or the classical case with no hopping. 
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The single particle spectral function can be denned as 



?M = --(S s l c * , ■ 1 p T7 
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where |gs) denotes the ground state, E is the energy of 
the ground state and r] the broadening constant. The 
calculation is done with a correction vector method. We 
keep up to 800 states and use a broadening constant of 
rj = 0.2. This rather large broadening saves computa- 
tional time and can be subsequently improved by per- 
forming a deconvolution of the spectral function^ 



where 1 is the unity matrix, J is the matrix of hopping 
amplitudes, and V = £j<5jj represents the matrix of dis- 
ordered on-site potentials. Inverting Eq. [7] yields the in- 
teracting local Green's function 

G lia (iu n ) = [fi a +ioj n -Ei-Y. lia (iuj n )-A ia .(iuj n )}~ 1 . (9) 

A set of new on-site hybridization functions &i a (iu> n ) can 
be extracted from the diagonal elements in Eq. [HI These 
constitute a set of new impurity problems to be solved 
by NRG, which closes the self-consistency loop. 



IV. RESULTS 



B. R-DMFT 

R-DMFT is the real-space extension of DMFT which 
incorporates site-dependency of the self-cncrgyi 6 ' 48 i 50 
Besides the description of the Mott-Hubbard metal- 
insulator transition and magnetic ordering it is capa- 
ble of treating spatial inhomogeneities such as disor- 
der. Each lattice site is mapped onto a single-impurity 
Anderson Hamiltonian Ha— 



la 
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where a\ a {cQ a ) and aj tT (co, T ) are fermionic creation and 
annihilation operators in the bath (on the impurity) and 
a represents the spin index. The parameters ei a and Vi a 
determine the hybridization function 
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which depends on the site index i and is determined self- 
consistently. 

For a given set of arbitrary hybridization functions, the 
solution of these effective quantum impurity problems 
is provided by the Numerical Renormalization Group 
(NRG) for T = and leads to a set of one- and 
two-particle on-site Green's functions, Gu (T (iu) n ) and 
Fu a (iuj n ) respectively. These determine the self-energy 
matrix in real-space representation^ 
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Due to the lattice Dyson equation, the interacting lattice 
Green's function is given by 



G a (iuj n ) 1 = G%(iu> n ) 1 - T, a {iui n ) 



(7) 



where w„ are Matsubara frequencies. The noninteracting 
Green's function G® a (iu) n ) in real-space representation is 
given by 



G°(iw n ) 1 = {ft a + iu n )l 



V. 
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A. Local density of states 

In order to describe the transition from dclocalized to 
localized states it is useful to characterize the spectral 
properties of the system by the local density of states 
(LDOS), which measures the local amplitude of the wave 
function at a given site i 



Pia{u) = lmG ila (u) . 
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The calculation of the on-site Green's function Gu a (uj) 
is based on the local spectrum for a particular disor- 
der realization {ei,E2, ■ ■ ■ ,£l}- In this work we focus 
on the paramagnetic solution and drop the spin index a 
for readability. To gain realization-independent informa- 
tion, arithmetic and geometric averaging over the spec- 
tral functions is performed 
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given N disorder realizations and L sites. In this work 
averages over 50 — 100 configurations were performed for 
each (U,D) parameter set. The geometric average is a 
good approximation of the typical value of the proba- 
bility distribution function for the LDOS 4^ Thus, the 
geometrically averaged spectral function is critical at the 
Anderson transition, i.e. it is finite in the delocalized 
regime and vanishes in the localized phase and there- 
fore can be interpreted as an order parameter.— i 44 ] 58 ; 59 
The arithmetic average remains finite in the Anderson- 
localized regime and corresponds to the global density of 
states (DOS). 

This classification is, however, only reliable in the ther- 
modynamic limit L — > oo. In finite systems a careful 
analysis of finite-size effects has to be performed. In 
a system of length L the spectrum is discrete and the 
non-interacting energy level spacing in the DOS scales 
as 1/L. This discrete level structure due to the finite 
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FIG. 1. pg(0) calculated by means of DMRG (upper panel) 
and R-DMFT (lower panel) . The Mott phase is located in the 
lower right area and the Anderson-Mott phase in the upper 
left area. They are separated by a region with finite typical 
LDOS, which indicates a derealization tendency. 



size can be smoothened by a broadening of each level to 
reconstruct the DOS in the thermodynamic limit. Un- 
fortunately, spectral broadening with a width 77, which is 
necessary due to the finite system size, limits the spec- 
tral resolution of our calculations. The interaction-driven 
metal-insulator transition can therefore only be detected 
when the gap exceeds 77. Similarly, the effect of the dis- 
order strength is underestimated. Disorder increases the 
energy level spacing in the LDOS, which leads to local- 
ization. Due to broadening this discretization is smeared 
out and Anderson localization is observed at a larger 
value of D than in the thermodynamic limit. A further 
finite-size effect is the competition between the system 
size L and the localization length £. When localization 
sets in with increasing disorder, the localization length is 
larger than the system size and even exponentially local- 
ized states contribute spectral weight at all lattice sites 
in finite-size systems. Accordingly, disorder-driven local- 
ization can only be identified for L ^> £. Due to these 
limitations, the goal of our work is to reveal the localiza- 
tion/delocalization trends in the system rather than to 
determine sharp phase boundaries. 

To investigate a possible emergence of a metallic phase, 
Pg(cj) is analyzed at the Fermi level u> = using R-DMFT 
with L = 128, 77 = 0.05 and 10-50 disorder configurations, 
and using DMRG with L = 64, 77 = 0.2 and 16 disorder 
configurations, where the results are additionally decon- 
volved. The results are presented in Fig. [TJ 
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FIG. 2. Comparison between geometrically (dashed line) and 
arithmetically averaged (solid line) LDOS for U = 3 calcu- 
lated within R-DMFT. With increasing disorder strength the 
geometrically averaged density of states gradually vanishes, 
starting from the band edges. 

Three different regimes are found in the system: The 
Andcrson-Mott insulating regime at large D, the Mott- 
insulator at large U and a derealization (metallic) ten- 
dency for intermediate disorder and interaction strength. 

(i) Anderson-Mott Insulator (AMI): As the spec- 
trum of the localized system consists of a dense distri- 
bution of poles^S the geometrically averaged LDOS van- 
ishes at all frequencies while the arithmetical average re- 
mains finite. The formation of AMI proceeds from the 
edges of the band towards the band center with increas- 
ing disorder strength (see Fig. [2]). As soon as the geo- 
metrically averaged LDOS p g {u) vanishes at the Fermi 
level oj = the system becomes fully localized. We first 
analyze our results along the D axis at U = 0. Exact cal- 
culations predict the non-interacting system to undergo 
an Anderson metal-insulator transition in the thermody- 
namic limit L — > 00 at any finite disorder strength^ 
However, we obtain clear signatures of localization at 
U = only for D > 2 (DMRG) or D > 3 (R-DMFT). 
This can be traced back to the finite size of the inves- 
tigated system. The localization length in this case is 
larger than the system size, such that the states remain 
quasi-extended and are only localized for very large sys- 
tem sizes j^i A careful analysis of the scaling behavior 
with system size was performed in R-DMFT for U = 
and D = 3 and is displayed in Fig. [3^- The decay of p g (0) 
with increasing system size confirms localization (see a 
detailed discussion of the finite-size scaling below). 

(ii) Mott Insulator (MI): Here for D = and U > 
the LDOS p(oj) vanishes at the Fermi level oj = due 
to the appearance of a Mott gap. In the deconvolved 
DMRG results p decreases monotonously with increas- 
ing interaction and a MI appears for U > 3. Within 
R-DMFT the Luttingcr theorem ensures pinning of the 
spectral weight at the Fermi levels Thus, p g {uj = 0) re- 
mains finite and constant up to U = 4 and the gap is first 
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FIG. 3. Finite-size scaling for the geometrically averaged 
LDOS at the Fermi-level calculated within R-DMFT. The re- 
sults are averaged over 10-50 configurations, (a) The decay 
with system size at D — 3 and U = indicates Anderson 
localization while for D = 3 and U — 3 (5) the vanishing 
size dependence denotes derealization, (b) Similarly, states 
which are Anderson-localized at 17 = become delocalized 
for (7 = 5. Finite interaction shifts p g (0) to higher values and 
reduces the dependence on the system size. 



formed at U > 5 (see Fig.HJ. DMRG as well as R-DMFT 
results for finite system size disagree with exact calcula- 
tions for a homogeneous system in thermodynamic limit, 
where a Mott insulator is predicted to appear at any fi- 
nite interaction strength^ In R-DMFT, as in every other 
DMFT extension, the finite metallic phase in ID is due to 
the approximation of a local self-energy, which becomes 
exact only in infinite dimensions. In DMRG this discrep- 
ancy in the detection of the Mott gap is a consequence 
of the unavoidable spectral broadening rj. Due to this 
broadening the Mott-gap is smeared out and the result- 
ing critical U represents an upper bound at which the 
Mott-insulator is formed. In DMRG, the broadening ef- 
fects are partially corrected by applying a deconvolution 
to the spectra. A careful analysis of the finite-size effects 
in the gap formation is presented in section ["IV Bl 

(hi) Derealization: For low to intermediate disor- 
der and interaction, the geometrically averaged spectral 



density p g {u}) is finite at the Fermi level w = 0. The non- 
black area in both phase diagrams in Fig. [1] indicates this 
derealization regime. Due to the deconvolution proce- 
dure within DMRG the spectral weights are lower than 
those obtained by R-DMFT. 

In the absence of disorder and interaction (U = D = 0) 
the system is metallic and satisfies Luttinger's theorem^ 
in agreement with our findings. Additionally, a der- 
ealization region is observed for intermediate interaction 
and disorder strengths. The latter, however, may in prin- 
ciple be due to overestimated spectral weights in finite 
systems for localization length £ > L and finite spectral 
broadening. 

To decide whether true derealization is observed, a 
careful finite-size scaling analysis of the R-DMFT results 
was performed for the intermediate U and D regimes. 
The scaling behavior of p g (0) characterizes the phase of 
the system: in the non-interacting case if the states are 
Anderson-localized and the localization length £ <C L, 
the geometrically averaged LDOS scales as exp[— L/£\ for 
periodic boundary conditions. If£ > L the decay of p g {0) 
with system size is algebraic^ In systems with finite 
interactions the functional dependence of p g (0) on system 
size is not known, however, the geometrically averaged 
spectral weight is still expected to decay to zero with 
increasing system size for localized states and to remain 
finite for delocalized ones in the limit L — > oo. 

In Fig. [3^ we observe the delocalizing effect of inter- 
actions. For D = 3 and U = the geometrically av- 
eraged LDOS decays with increasing system size, which 
indicates the Anderson localized phase. Thus, the final 
spectral weight visible in the phase diagram in Fig. Q] 
(lower panel) for this parameter set is a finite-size effect 
for the L = 128 lattice. Upon increasing the interaction 
strength to U = 3 and U = 5 no clear localization signa- 
ture could be observed anymore within accessible system 
sizes. 

In Fig. [HJd the effects of increasing disorder strength 
and interactions on the localization properties of the sys- 
tem are compared. A geometric average which decays 
with increasing system size L for U = and D = 2 and 3 
indicates the Anderson localized phase. Again, the der- 
ealization trend observed in Fig. [T] (lower panel) for this 
parameter set turns out to be a finite-size effect for the 
L = 128 lattice. The redistribution of spectral weight due 
to increased disorder strength leads to a well pronounced 
shift of p g (0) towards lower values for D = 3 compared 
to D = 2 in the non-interacting case. Strong interaction 
U = 5 shifts the geometric average upwards and the de- 
pendence on the system sizes studied here vanishes which 
indicates derealization. 



B. Charge gap 

In the homogeneous ID Hubbard model at half-filling 
a metallic phase exists only at U = 0. At any finite in- 
teraction strength a charge gap in the density of states 



0.01 0.02 0.03 0.04 

1/L 



1 



2 
D 



3 4 



FIG. 4. Charge gap versus inverse system size for U — 3 
calculated within DMRG. For disorder strength D > 1.5 the 
disappearance of the gap in the thermodynamic limit 1/L — > 
indicates the transition from the gapped Mott insulating to a 
gapless delocalized phase. 



FIG. 5. Scaling of the gap normalized to the non-disordered 
values Go in the thermodynamic limit, calculated within 
DMRG. The gap effectively vanishes for disorder strength 
D ~ {7/2. The error bars correspond to the statistical er- 
rors of the fitting routine. 



is predicted to appear and the system becomes Mott in- 
sulating. Thus, starting in the MI phase, the additional 
influence of disorder can be detected via a vanishing gap. 
The charge gap is calculated by means of DMRG as 
G = [E(N + 1) + E(N - 1) - 2E(N)}, where E(N) is 
the ground state energy for N particles. It is important 
to note that the computation of the energy E(N) is not 
based on the spectrum and is therefore not affected by 
the artificial broadening ?/. Thus, the calculated charge 
gap is exact for a given system size. We use system sizes 
up to L — 128 sites and average over 16 disorder realiza- 
tions. 

The finite-size scaling for the charge gap at U = 3 
and various disorder strengths is presented in Fig. |4j In 
the thermodynamic limit 1 / L — > the charge gap closes 
with increasing disorder and vanishes at D w 1.5. The 
scaling of the charge gap with system size was fitted by 
the function G(L) = Goo + fj- This finite-size analysis 
confirms the derealization in Fig. [1] (upper panel) for 
[7 = 3 and D > 1.5. In contrast, the finite spectral 
weight at smaller disorder strength was found to be a 
finite-size effect. 

The extrapolated gaps for various disorder and inter- 
action strengths normalized with respect to the extrapo- 
lated gaps of the corresponding homogeneous systems are 
presented in Fig. [S]for different interaction strength. The 
indicated error bars are the asymptotic standard errors of 
the fitting routine. For intermediate interaction strength 
the critical disorder needed to destroy the gap was found 
to be ps U/2. This tendency becomes less clear for U < 2 
as the limited resolution of the exponentially small gap 
complicates our analysis. The remaining gapless region 
covers the delocalized as well as the Anderson-Mott lo- 
calized regime. For a quantitative analysis of this region 
we calculate the inverse participation ratio as described 
in the following section. 
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FIG. 6. (Color online) Finite-size scaling of the IPR, calcu- 
lated within R-DMFT. The U = 0, D = 3 data are fitted by 
1(0) = Zoo(w) coth[L/£(ti;)] which corresponds to the charac- 
teristic scaling in the localized regime. The decay of 1(0) for 
D = 3 and V — 3 (5) with increasing system size indicates 
derealization. 



C. Inverse participation ratio 

The localization properties of the system can be quan- 
tified by means of the inverse participation ratio I(uj) 
(IPR). This observable corresponds to the inverse of 
the number of sites over which a state is extended. In 
the non-interacting homogeneous case it is defined as 
I(oj) = Y^f l ? /'i( w )| 4 i where ip(to) is a single particle wave 
function. If a state at frequency u is exponentially lo- 
calized with corresponding localization length £(w) then 
I(u>) scales as ioo(w) coth[L/£(uj)] with the system size 
for periodic boundary conditions and lim£_ i . 0O /(cl>) = 
-Too (<*>) = 1/C( W )- In the case of a purely delocalized 
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state, i.e. for a wave function homogeneously extended 
over the whole system, the IPR vanishes as 1/L with 
increasing system size. 



In an interacting system the IPR is defined a; 
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(13) 



Similarly to the non-interacting case, in the thermody- 
namic limit the delocalized states lead to a vanishing 
IPR, while in the localized system the IPR remains fi- 
nite for L —> oo. Thus, the scaling behavior with system 
size gives insight into localization or derealization driven 
by disorder or interaction. In Fig. [5] the IPR at us = 
is plotted as a function of system size. As expected, for 
U = and D = 3 the system is Anderson localized and 
the IPR saturates for L —> oo. Moreover, fitting the data 
by 1(0) = loo coth[L/£(w)] shows that the scaling be- 
havior also holds for open boundary conditions at larger 
system sizes. For finite interaction strength U = 3 and 
U = 5, however, 1(0) decays with system size (note the 
differing scales on the y axis). The vanishing IPR high- 
lights the dclocalizing effect of the repulsive interaction. 
These results confirm our findings in Fig. |3l 



V. SUMMARY 

We have studied the occurrence of derealization, 
Anderson- and Mott-insulating phases for interacting dis- 
ordered spinful fcrmions on a ID lattice at commensurate 
filling. The phases were characterized by means of the 
local density of states, the scaling of the inverse partic- 
ipation ratio and the charge gap. Both numerical tech- 
niques used here, DMRG and R-DMFT, agree in their 
predictions for parameter regimes well within the Ander- 
son and Mott phases. However, R-DMFT was found to 
be more affected by finite-size effects than DMRG, which 
is reflected in a larger derealization trend for fixed lat- 
tice size. As our main result, we present the first non- 
perturbative calculations indicating a delocalized phase 
of spinful fcrmions in ID due to the interplay of disor- 
der and interaction of intermediate strength, similar to 
the bosonic case. This trend, observed in DMRG as well 
as R-DMFT, persisted in finite-size scaling within both 
methods. 
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